Hole dynamics in spin and orbital ordered vanadium perovskites 
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Hole dynamics in spin and orbital ordered vanadates with perovskite structure is investigated. A 
mobile hole coupled to the spin excitation (magnon) in the spin G-type and orbital C-type (SG/OC) 
ordered phase, and that to the orbital excitation (orbiton) in the spin C-type and orbital G-type 
(SC/OG) one are formulated on an equal footing. The observed fragile character of the (SG/OC) 
order is attributed to the orbiton softening caused by a reduction of the staggered magnetic order 
parameter. It is proposed that the qualitatively different hole dynamics in the two spin-orbital 
ordered phases in vanadates can be probed by the optical spectra. 



PACS numbers: 71.30.+h, 71.10.-w, 75.30.-m, 78.30.-j 



One of the central issues in the correlated oxides 
is doped Mott insulator [1,2]. The prototypical com- 
pound is the high Tc superconducting cuprates (HTSC) 
where the parent Mott insulating state with the anti- 
ferromagnetic (AFM) long-range order is understood by 
the single-band electronic model with strong Coulomb 
interactions. Apart from the half filling, doped holes 
are strongly renormalized by interaction with spin ex- 
citations. 

There is another class of Mott insulator where the or- 
bital degeneracy of the transition-metal ion remains. Mo- 
bile carriers doped into this class of insulator couple to 
the spin and orbital arrays. Vanadates with the per- 
ovskite structure Ri-xAxVO^ {R: a trivalent rare-earth 
ion. A: a divalent alkaline-earth ion) are the likeliest can- 
didate. A number of researches have been done in con- 
nection with the spin-orbital orders in insulating i?V03 
[3-12]. In a V^"*" ion, two electrons occupy the t2g or- 
bitals with the S — 1 high-spin state. Two kinds of spin- 
orbital orders are found in the ground state [3,5]: the 
G-type (three-dimensional (3D) staggered) AFM order 
associated with the C-type (rod type) orbital one termed 
the (SG/OC) phase in YVO3, and the C-type AFM or- 
der with the G-type orbital one termed the (SC/OG) 
phase in LaVOa. Both the orders have the commonly 
occupied d^y orbital at each V site, and the alternately 
occupied dyz/dzx orbitals with the AFM spin alignment 
in the ab plane. A difference between the two is seen 
along the c axis; the staggered spin (orbital) order associ- 
ated with the uniform orbital (spin) one in the (SG/OC) 
((SC/OG)) phase. 

Apart from the well studied i?V03, little is known 
about nature of mobile holes in the spin-orbital ordered 
Ri-xAx^O^. A dynamics of a d^y hole is only to be 
expected from the consequence in a well-known doped 
two-dimensional (2D) quantum AFM and is common in 
the (SG/OC) and (SC/OG) phases. A clear contrast 
is anticipated in a hole motion along the c axis; in the 
(SG/OC) and (SC/OG) phases, a mobile hole is strongly 
coupled to the spin and orbital excitations, i.e. magnon 
and orbital, respectively. Orbiton, standing for a quan- 



tized object of the collective excitation in orbital order, is 
observed in the Raman scattering of LaMnOs [13] and is 
also proposed in RVO3 [12]. Thus, the perovskite vana- 
dates are idea for study of the quantum hole dynam- 
ics in a spin-orbital ordered states. One of the striking 
difference between the observed two phases is the sta- 
bility; the (SG/OC) phase is rather fragile against the 
hole doping [4,14]. the (SG/OC) phase realized in YVO3 
rapidly disappears by doping of holes about x < 0.05 
in Yi-^Ca^^VOa. This phase is changed into (SC/OG) 
which is stable, at least, up to a; '--^ 0.4. The (SC/OG) 
phase in a series of Lai_a;Sra;V03 also survives until 
X = 0.178 [5]. This fragile nature of (SG/OC) is not 
explained within the doped 3D AFM with S = 1. 

We present, in this letter, a theory of doped perovskite 
vanadates with spin-orbital orders. The hole dynam- 
ics associated with magnon and orbiton in the (SG/OC) 
and (SC/OG) phases are formulated based on the self- 
consistent Born approximation (SCBA). The different 
quasi-particle (QP) properties between the two phases 
and the observed barely stable (SG/OC) phase are at- 
tributed to the softening of the orbital excitation. 

Our starting point is the t — J type model with the 
t2g orbital degree of freedom, Ti = Tit + Tij + Hjt, de- 
rived from the generalized Hubbard model in the large 
on-site interaction limit [15]. The exchange term between 
nearest neighbor (NN) spins and orbitals has the form 
= n^A^ + 'H-'E + n2T, + H2T^ where ^^2, ^E, ^Ti, 
and indicate the intermediate {t2gY configurations 
of the exchange processes. The leading term involves the 
high-spin state explicitly given as 
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where Si is the spin operator with magnitude S = 1. 
JiA2 {— (U' — I) = J) is the exchange interaction with 
the on-site inter-orbital Coulomb interaction U' , the ex- 
change interaction / and the transfer intensity t between 
NN V ions. The other exchange interactions are ex- 
pressed by J and rj = I /U' . The orbital sector in Tiij^ 
involving A', i?', and C', where /(= x,y,z) indicates 
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a direction of a bond, is expressed by the local orbital 
operators at sites i and j. These describe the diagonal 
and off-diagonal components of the electric-quadrupole 
moments, OiE-f (7 = u.v) and OiT27 (7 = z), re- 
spectively, and the orbital angular moments OiT^ (7 = 
x,y,z). There is no continuous symmetry in the orbital 
sector of Hj. The t term is given as 



(2) 



The operator Cj-yo- annihilates a t2g electron at site i with 
orbital 7 = {yx, zx,xy) and spin a = (T,i), and is de- 
fined in the Hilbert space where the (t2g)"^ configuration 
is excluded at each site. In the ideal perovskite crys- 
tal, the hopping between the different orbitals is prohib- 
ited, and one of ijj^s is zero and other two arc equal, 
e-g- *M+z = and = tl^^_^_^ = t. In addition to 

the t and J terms, roles of the Jahn- Teller (JT) inter- 
action described by TijT = 9^1^ QiE-yOiE-f are stressed 
in the insulating RVO3 [9]. QiE-y's (7 = u,v) are the 
normal modes of an Oe octahedron with symmetry Egj 
and are treated as classical valuables. The detail expres- 
sion of the above Hamiltonian and studies of the orbital 
order/excitation are presented in Refs. [12] and [15]. Sim- 
ilar spin-orbital models for insulating KVO3 are derived 
by other authors [7-10,16,17]. 

Dynamics of holes doped into the spin-orbital ordered 
vanadates are formulated by SCBA which has been ap- 
plied to HTSC [18,19] and the colossal magnetoresistive 
(CMR) manganites [20-22] . The electron operator is ex- 
pressed in the slave-fermion representation in the large 
limit of S. We assume that the orbital excitation between 
the dyz and dzx orbitals couples effectively to the doped 
holes, since (i) the excitation between d^y and dyz /dzx 
dose not propagate coherently (local modes) [12], and (ii) 
the dxy band is separated from the dyz /dzx ones, i.e. the 
excitation energy between dxy and dyz /dzx is higher [6]. 
The electron operator is rewritten as. 



^7={xy) filler + ^'1=(zx,yz)fiSi„t^^ 
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We introduce two fermionic operators {fi, fi} for charge, 
two bosonic operators {sic^jSic^} for spin, and bosonic 
one {U^} for orbital with the local constraints fjfi + 
Ea«Lsi<^ = 1' fifi + Ea4crSi<^ = 1 and E^sL-S':- = 
J2-y=iyz.zx) A-y^ii- The Operator Ps projects onto the '^Ti 
high-spin state in the (^23)^ configuration. Then, Tit is 
decoupled into the dyz/ dzx orbital part and the dxy one 
as Ht = Uyzizx + Uxy with Uyzizx = ^^^^ + H^'^^'l Thc 
main term is Ti^'^^ describing a motion of the dyz/ dzx hole 
along the c axis of our present interest. In the large limit 
of 5*, we remain the down (up) boson operators in the up 
(down) spin/orbital sublattice. Ti'^'^-' has the form 
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FIG. 1. Contour plot of spectral function and DOS of elec- 
tron in (SG/OC) ((a) and (b)), and those (SC/OG) ((c) 
and (d)) for a; = 0.1. Broken lines in (b) and (d) are 
for A{k = (0,0,^), w) multiplied by 1/5. Parameters are 
J/t = 0.2, rj = 0.125 and qQe/J = 1.25. The origin of the 
vertical axis indicates the Fermi energy. The inset shows the 
J dependence of the QP weight at A; = (0, 0, ^). 



where = Sj^ for the (SG/OC) phase and uj: ~ for 
(SC/OG) with g(k{) — cosafc;. Explicitly, the mobile 
hole along the c axis is coupled to magnon (orbiton) in 
the (SG/OC) ((SC/OG)) phase. Ti^"*) andW^^, are com- 
mon in the two phases. Ti^"^-* is for the dyz/ dzx hole in 
the ah plane, and fomrs f-fl stit - The simulta- 

' ■' y k' 1 k'-k+q 

neous magnon-orbiton emissions/absorptions is expected 
to suppress the QP weight and dispersion. Tixy is for 
the dxy hole in the ab plane, Hxy = J2k,k'{9{kx) + 

9{kv)}f-^f\,^k-k' ~^ H.C.. Similar physical consequences 
arc only to be expected from this term to those in a 
mobile dx^-y^ hole in the 2D quantum AFM studied in 
HTSC [18]. Here we focus on H^'^^ characterizing differ- 
ent hole dynamics in the (SG/OC) and (SC/OG) phases. 
In the same scheme, Tij is formulated by the spin wave 
expansion [12]. In the magnon (orbiton) part of Tij 
written by Si^ (ti-y), the orbital (spin) operators are re- 
placed by the static correlation functions {OiT-yOjT'Y) 
{{Si ■ Sj)), and the simultaneous spin-orbital excitations 
arc neglected, same as those in Tli mentioned above. The 
lattice effects arc limited to be static ones given by H-jt, 
although the coupling with phonon may modify the QP 
band width and the orbiton dispersion around the cross- 
ing point with phonon in the small J limit. 

The electron spectral function A{q,ui) is calculated 
from the fermion itor G{k, llj) where the boson op- 

erators are replaced by the saddle-point solutions. The 
self-energy for G{k, ui) is given by 

^{k,u;) = j dp^G{k- q,uj-u) 



X \gu{k,q)\^Duo{q,i^), 



(5) 



where u = s for (SG/OC) and u = t for (SC/OG). The 
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FIG. 2. Contour plots of the spectral functions at a; = 0.1 
and the doping dependence of DOS for orbiton ((a) and (b)) 
and those for magnon ((c) and (d)) in (SG/OC). White bro- 
ken lines in (a) and (c) indicate dispersions at a; = 0. The 
lower branch of the orbiton in (a) corresponds to that along 
(0,0, ^) - (0,0, ^) - (^,0, 1) in the Brillouine zone for the 
cubic unit cell. Parameters are the same with those in Fig. 1. 



bare magnon/orbiton propagators Du{){q,v) arc used in 
Eq. (5) instead of the full propagators Du{q, v)- gu{k, q) 
is the coupling constant including the coefficients in 
the Bogoliubov transformation. A{k,uj) along (0,0,0) — 
(0, 0, ^), with a being the cubic perovskite unit cell, and 
the density of states (DOS) N{ui) are obtained in the self- 
consistent calculation (Fig. 1). A reasonable parameter 
set as J/t = 0.2, rj = 0.125 and gQs/J = 1 ~ 1.25 with 
J being around SOmeV is chosen based on the previous 
evaluations [9,12]. In both the two phases, A{k,uj) shows 
the broad incoherent continuum and the well-separated 
QP peak near the Fermi level (FL). The dispersive QP 
peak has the band width D of the order of J and has the 
highest energy at the momentum (0,0,^). The clear 
differences between the two phases are seen in the QP 
spectral weight a(k) and D; in the (SC/OG) phase, a 
large portion of the spectra are concentrated on the QP 
peak with the larger width D, and a{k) is rapidly sat- 
urated with increasing J (the inset of Fig. 1). This is 
in contrast to the case of (SG/OC) where the spectra 
are dominated by the incoherent components. These are 
ascribed to the larger excitation energy of orbiton with 
the excitation gap about 2J ~ 3J (as shown later in 
Fig. 3(b)). We have checked and confirmed these fea- 
tures of a{k) and D by the variational-typc calculations 
[23]. It is worth noting that the hole dynamics associ- 
ated with the t2g orbitons in (SC/OG) is in qualitatively 
contrast to that with the Cg orbital for the CMR man- 
ganites [24,20]; A{k, co) has a large weight around the free 
eg bands rather than the QP part. This is because of the 
non-zero hopping matrices between the NN different Cg 
orbitals unlike the t2g ones. 

Now we turn to the spin and orbital sectors. The renor- 
malized magnon/orbiton propagators are given by utiliz- 
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FIG. 3. (a) Doping dependences of reductions of the spin 
and orbital order parameters, AS and AT, in (SG/OC). The 
inset shows those in (SC/OG) with J ft = 0.2. (b) Doping de- 
pendence of the minimum orbiton energy uJow- The filled and 
open symbols in (b) are for qQe/J =1.25 and 1, respectively. 



ing the calculated G{k,uj). The bubble-type diagram is 
the lowest order of the self-energy, 



n«(g,a;) = 



2t:N 



X G{k + q,uj + u)G{k, v). 



(6) 



The spectral functions S{q,uj) = (— l/7r)ImDs(<f, tj) and 
T{q,uj) = (-l/7r)ImA(5,w), and DOS in (SG/OC) are 
shown in Fig. 2. In the undoped case (white broken 
lines), the calculated S{q,u)) well reproduces the exper- 
imentally observed spin wave dispersion in YVO3 [11]. 
The flat dispersion of T{q,uj) in the ab plane originates 
from no coherent propagation of the dy-^/d^x orbital ex- 
citations in this plane [12]. The excitation gap in T{q, u>) 
at X = comes from the discontinuous symmetry of 
the orbital sector in Tij and the JT interaction, and is 
much smaller than that in (SC/OG). Namely, the frag- 
ile character is inherent in the C-type orbital order. As 
expected from distortions of the staggered spin/orbital 
array by mobile holes, softening and broadening are ob- 
served for Siq,u;) in (SG/OC) (Figs. 2(c) and 2(d)) and 
for T{q,uj) in (SC/OG) (not shown). Unexpectedly, the 
band width of T(g, w) in (SG/OC) becomes broader with 
doping (Figs. 2(a) and 2(b)), in spite that the orbital dis- 
tortion due to the hole motion is not expected along the c 
axis. The minimum energy Wqw -^(^' ^^^^ phase 
is located at the point F in the lowest branch which cor- 
responds to the point (0,0, ^) in the original Brillouine 
zone for the cubic perovskite unit cell. cOq^ decreases 
with doping and finally touches the zero energy at small 
number of x. That is, the (SG/OC) phase is unstable 
against the small doping and is changed into the G-type 
orbital phase. This orbiton softening in (SG/OC) origi- 
nates from reduction of the staggered spin order parame- 
ter, termed AS'(> 0) (see Fig. 3). The remarkable change 
by doping is only observed in AS of the (SG/OC) phase. 
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However, this reduction itself is not sufficient for instabil- 
ity of the AFM ordering, i.e. AS < S{= 1). Instead, the 
effective exchange interaction for the staggered (uniform) 
orbital ahgnment in Tij, given as J'^p = JiA2 (2+(<S'i-5^)) 
{Jf" = Jr(l - {Si ■ Sj)) (r =^E, 2Ti, 2T2)), rapidly 
increases (decreases) with increasing AS. Finally, the 
orbiton energy becomes zero at the critical x which de- 
pends on the exchange and JT interaction parameters. 
In short, the instability of the (SG/OC) phase with dop- 
ing, observed in Yi-^Ca^VOs [14], is caused by the or- 
biton softening due to increasing of AS through the spin- 
orbital coupling. 

The qualitatively different hole dynamics in the two 
phases can be probed by the optical experiments. The 
2;-component of the optical conductivity spectra (Tzz(a;) 
is formulated in the slave- fcrmion scheme with the 1/S 
expansion. The lowest order diagram of (Tzz{ijJ) is given 

by <Jzz{^^) = — ^"c2^ lmK{co) with 
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FIG. 4. Regular parts of the optical conductivity spectra 
for (SG/OC) (a), and those for (SC/OG) (b). The inset shows 
the integrated spectral weights for the (SG/OC) phase (bro- 
ken lines) and those for the (SC/OG) ones (bold lines). Pa- 
rameter values are the same with those in Fig. 1. 



merical calculation has been performed by the supercom- 
puters in IMR, Tohoku Univ., and ISSP, Univ. of Tokyo. 



X Duo{k2 — ki,u; - u!i + W2)G{k2,oJ2), 



(7) 



for u = s {u = t) in the (SG/OC) ((SC/OG)) phase, 
indicating the particle-hole pair creation associated with 
magnon/orbiton. Here we introduce the results of the 
regular part of azzi^jj) at finite frequency, rather than the 
Drude part (Fig. 4). With increasing holes, the sharp 
peak structure growing up around a; = J 2J in the 
(SC/OG) phase, reflecting the large QP weight in A{k, uj) 
(Fig. 1), in contrast to the broad incoherent spectra in 
(SG/OC). A shoulder structure is clearly seen in the inte- 
grated spectral weight M{uj) = azz{^^'')dcj' (the 
inset of Fig. 4) which would be an evidence characterizing 
the hole dynamics in the orbital ordered insulators. 

In summary, we present a theory of doped vanadium 
perovskites with spin and t2g orbital orders. The doped 
(SG/OC) and (SC/OG) phases are characterized by the 
hole motion along the c axis scattered by magnon and 
orbiton, respectively. The QP dynamics associated with 
the t2g orbital excitation show remarkable difference from 
that with magnon as well as the Cg orbital excitation. In- 
stability of the (SG/OC) phase into (SC/OG) observed 
in Yi_a;Caa;V03 is attributed to the orbiton softening 
through a reduction of the staggered spin order param- 
eter. The angular resolved photoemission spectroscopy 
and optical experiments for the QP, as well as the inelas- 
tic neutron scattering for magnon and orbiton [12], are 
required to test directly the present theoretical results. 
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